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Abstract 

The mass and decay width of the © + (1540) with isospin are investigated in a constituent quark 
model comprising uudds quarks. The resonance state for the + is identified as a stable solution 
in correlated basis calculations. With the use of a one-gluon exchange quark-quark interaction, 
the mass is calculated to be larger than 2 GeV, increasing in order of the spin-parity, | , | and 
^ + (| + ), and only the | state has a small width to the nK* + decay. If the calculated mass is 
shifted to around 100 MeV above the N + K threshold, the 6+(1540) is possibly ^ + (| + ) or |~, 
though in the latter case it cannot decay to the nK + channel. In addition it is conjectured that 
other pentaquark state with different spin-parity exists below the O + (1540). The structure of 
the + is discussed through the densities and two-particle correlation functions of the quarks and 
through the wave function decomposition to a baryon-meson model and a diquark-pair model. 
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I. INTRODUCTION 



Quark model succeeds to classify the hundreds of baryons and mesons. A constituent 
quark model well reproduces the mass spectra of these hadrons using a colored interaction of 
one-gluon exchange (OGE) type together with a phenomeno logical confinement potential 0, 
13]. There are of course those hadrons which do not fit the prediction of the constituent quark 
model, e.g., iV(1440) (Roper resonance) and A(1405). Recently, motivated by a theoretical 
prediction for an exotic baryon [3], some experimental groups have searched for it and 
reported the observation of G + (1540) which has baryon number +1, strangeness +1 and a 
decay width of less than 25 MeV [3, 13, El- Its minimal quark content is uudds, so it is often 
called a pentaquark baryon. A reason why the + attracts much interest is that it is truly 
exotic, its mass is about 100 MeV above the N+K threshold, and its decay width appears 
to be very small nevertheless. The spin and parity of the + has not yet been known, nor 
even the existence of the 9 + (1540) seems to be established experimentally 0,13]. 

The recent controversial status of the O + (1540) warrants a careful study of this system. 
From a theoretical side a calculation of the pentaquark mass and its decay width is useful 
and important. It was speculated that the + is a bound state of two highly correlated 
ud pairs (diquark) and s [3] or it is a coupled system of a color antisymmetric ud diquark 
with a uds triquark [13]. In these models the spin and parity of the G + was assumed to be 
J P = \ to account for its small decay width. 

Some dynamical calculations of the mass of the + for different J p states have been 
performed in constituent quark models 11, 13, 13, 14 1 . The quark model usually predicts a 
smaller mass for the negative parity G + than for the positive parity G + . In fact the ground 
state of the + is ~ , but it was argued that the mass difference between the \ and | + 
states may become small or even be reversed for a certain type of flavor-spin dependent 
interactions acting between the quarks [13]. See a review article [HI for some theoretical 
attempts and problems which are discussed from a broad perspective. 

All the calculations performed so far in the constituent quark model were based on a 
variational principle which can be applied to a bound state problem. Without giving a link 
between such a bound state calculation and a decay width, it is hard to quantify the width 
of the + . Though there are some attempts to predict the decay width of the + in the 
constituent quark model jl2l 17, 13], either the unbound nature of the G + is not taken into 



account or the assumption of its structure is too simple. Moreover, the correlated motion 
among the five quarks is usually truncated in the orbitaL spin, isospin and color degrees 
of freedom or in some of them apart from a calculation |l3[. For example, a treatment of 
hidden color components is not quite clear. Though the colors of the quarks are uniquely 
coupled to a color singlet state in a baryon and a meson, we have several possibilities to 
get a color singlet pentaquark. Therefore it is desirable to carefully examine the role of 



the hidden color states in a variational calculation. This problem is avoided in Ref. [ll[ by 
averaging a confinement potential. 

We will instead show that it is possible to predict the mass of the pentaquark even 
though the hidden color components are included in a calculation. A very recent quark 
model calculation with a scattering boundary condition [14| is perhaps the only one that 
has undertaken a prediction of the width of the G + . It is claimed there that a sharp resonance 
with \ appears above 2 GeV, which is in disagreement with the other predictions 13, 13] • 
In any case, the appearance of such a narrow resonance is hard to understand unless the + 
has negligibly small NK component. It is fair to say that the quark model study has not 
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so far given a settled prediction for the decay width of the + though it seems to predict 
more or less a similar mass for it. 

The purpose of this study is to predict the mass and the decay width of the 9 + in a 
few-body approach based on the quark model. We assume that the + is a five-particle 
system but not a heptaquark system like KtcN as discussed in 2(| . We calculate 



the mass and the width for some J p states assuming that the system has isospin and is 
confined in a spatially small region. We pay a special attention to the unbound nature of 
the sought solution and perform a dynamical calculation which allows for not only spatial 
correlations among the quarks but also all possible spin-isospin-color configurations. The 
merit of employing the quark model is that it exploits the symmetry of the system and can 
test the validity of some models through the dynamical calculation. In this way, we can 
learn the extent to which the diquark model |9| or the diquark-triquark model [lCj is sound. 
Also our model can include the effect of NK and NK* channels on the + , which is vital 
to predict the width. 

We use a correlated Gaussian basis and a complete set of the spin-isospin-color configura- 
tions. To predict the mass and decay width of the + , we use a real stabilization method 
which utilizes square-integrable basis functions to localize the resonance. We use a stochastic 
variational method 0, (SVM) to set up the basis set which describes resonance states. 
Our emphasis here is not to include continuum states explicitly but to localize the resonance 
in a simple way, which makes it possible to calculate the reduced width. We believe that 
to evaluate the decay width is still a worthwhile study, considering the present controversial 
situation as mentioned above. 

The composition of the paper is as follows. Our formalism is explained in sect. |HI We 
use an OGE potential as the interaction between the quarks, as defined in sect. Ill Al Some 
details are given for the basis function in a particular representation in sect. Ill Bl and for its 
transformation to other basis functions in sect. Ill 01 The application of the real stabilization 
method to the description of a resonance is presented in sect. IIHI Here we illustrate an 
emergence of a unique solution (resonance) in the bound-state-looking calculations. The 
results of calculation are presented in sect. IIVI The mass spectrum of the + is given in 
sect. II V Al the decay width is discussed in sect. II V Bl and the structure of the 6 + is analyzed 
in sect. II V CI A summary is given in sect. |3 



II. FORMALISM 
A. Hamiltonian 

A Hamiltonian for the + reads as 

5 5 2 

^ = E^ + E|^-- T c. m . + E^ ) (!) 

i=l j=l zr 'h i<j 

where is the mass of the ith quark and Vy is the interaction potential between the 
quarks. The so-called natural units are used, so a length has a dimension of inverse energy 
(1 fm= I g^3 MeV" 1 ). The kinetic energy of the quark is given in a nonrelativistic form and 
the kinetic energy of the center of mass motion, T c m , is subtracted from the Hamiltonian. 
The eigenvalue of the Hamiltonian corresponds to the intrinsic mass of the + . 

The quark-quark potential used in this paper consists of a variant of OGE potential, 
a phenomenological confinement potential and a zero-point energy term. It is taken from 
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Both potentials contain the color Coulomb and color magnetic terms. 

Af are the color SU(3) generators (Gell-Mann matrices) for the ith quark. In this paper 
an SU(3) representation is labeled by Elliott's convention |24| T = (A/i) whose dimension 
is d(Xfi) = |(A + l)(/i + l)(A+/i + 2). A quark carries a color of (10) representation and an 
antiquark carries that of (01) representation. The color state of two quarks is either (20) 
(color symmetric) or (01) (color antisymmetric). Corresponding to these representations, 
the matrix of (Af • Af ) becomes 



(Af ■ Af ) 




(4) 



For the case of a quark-antiquark pair the color state is either (11) or (00), and the matrix 
of (Af • Af ) corresponding to these states is 



(A 



A?) 




(5) 



In cases where the two quarks are in the color state (20) or the quark-antiquark pair is in 
(11), (Af ■ Af ) gives a sign opposite to the physical cases which appear in a baryon and a 
meson, so those color states play a role of deconfming the quarks at large distances. We will 
not exclude such color states, however, from the beginning as they may be important in the 
@ + comprising four quarks and one antiquark in a spatially confined region. 

For a color singlet system, the zero-point energy term has an expectation value 
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(6) 



where the subgroup label of the color singlet state is denoted 000, and where n is the total 
number of quarks and antiquarks contained in the system and it determines the expectation 
value independently of the way of constructing the color state. Thus the zero-point energies 
of the AL1 potential are — 2A, — 3A and — 5A, for mesons, baryons and the + , respectively. 
In the case of the TS potential, v is varied depending on n: v = Vq for the meson, v = 3V 
for the baryon, and vq = 5Vq for the + . To change the parameter of the zero-point energy 
results in simply shifting the mass but never alters the eigenfunction of the Hamiltonian. 
This property will be used in sect. HVBI to estimate the decay width of the + . 
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We list in Table [Q the masses of some mesons and baryons calculated using the 
SVM jiil . 23[. The parameters of the potential are also given in the table. It is found 
that both potentials give the results which agree reasonably well with the observed masses 
of K, K*,N(939) and A(1232). The N + K threshold is 1486 and 1451 MeV for the ALl 
and TS potentials, respectively The masses of the strange baryons predicted with the TS 
potential are considerably large compared to the experiment, which is probably due to the 
fact that the mass of the strange quark is taken to be large. The masses of some baryons 
such as A^(1440) and A(1405) are found to be too large. A flavor-dependent potential or 
some multi-quark configurations may be needed to reproduce these masses ji^, |26|. The 
root mean square (rms) radius of the quark distribution is also calculated and listed in the 
table. The calculated size of A^(939) is apparently too small, so we must make allowance 
for this underestimation in setting a channel radius which is needed to calculate the decay 
width of the 9+ 



B. Basis function 

Our Hamiltonian commutes with the total orbital angular momentum, the total spin and 
the total isospin. Thus the + state is specified with their quantum numbers L, S and T. 
By letting J and P denote the total angular momentum and parity of the B + , we may write 



TABLE I: The masses M and rms radii a/ (r 2 ) of mesons and baryons. L is the total orbital angular 
momentum assumed in the calculation. The parameters for the ALl potential 0] are m U( i = 315 
MeV, m s = 577 MeV, k = 0.5069, A = 0.1653 GeV 2 , A = 0.8321 GeV, k' = 1.8609, 5 = 0.2204, and 
A = 1.6553 GeV 8-1 , while those for the TS potential [13] are m ud = 313 MeV, m s = 680 MeV, 
05 = 1.72, A 9 = 3 far 1 , o C onf = 172.4 MeVfar 1 , V = -345.5 MeV, and ^' = -742.7 MeV. 



ALl TS 



Particle L M [MeV] y/lr 1 ) [fm] M [MeV] yflff [fm] 



vr(138) 





138 


0.30 


106 


0.32 


p(770) 





769 


0.46 


660 


0.45 


K(496) 





491 


0.31 


514 


0.33 


K*(892) 





903 


0.42 


813 


0.41 


AT(939) 





995 


0.49 


937 


0.50 


A(1440) 





1722 


0.75 


1755 


0.73 


A(1232) 





1307 


0.58 


1229 


0.56 


A(1116) 





1148 


0.47 


1266 


0.49 


A(1405) 


1 


1522 


0.58 


1736 


0.60 


£(1193) 





1229 


0.49 


1381 


0.51 


5(1318) 





1349 


0.44 


1659 


0.49 


0(1672) 





1675 


0.48 


2084 


0.50 
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its wave function as an antisymmetrized product of the orbital, spin, isospin and color parts: 



JM 



A [^ orbital) 4 spin) ]jM«rv 



(isospin) J (color) 
" '(00)000 



(7) 



where A is an antisymmetrizer of the four quarks (uudd) labeled 1, 2, 3, and 4 (s is labeled 
5) normalized to A 2 = A, the square bracket [ , 0^ >rblt V's ] JM stands for the angular 
momentum coupling. The G + masses with different J values for a given set of L and S 
values are degenerate in our model. 

Now we specify each part of Eq. (JJJ) in detail. It is vital to allow for various types of 
correlation among the quarks. For the orbital part, we use an explicitly correlated Gaussian 
basis in a global representation (22!. I231 127^: 



(orbital) 
LMt. 



") LMl (A,u,x) = exp^--xAx^y LMt (ux) 

= exp j-i AijXi ■ x-^j y LML (Y^UiX^J 



(8) 



where x (x±, X2, x$, X4) stands for a set of the intrinsic Jacobi coordinates other than the 
center of mass coordinate of the system and it is defined from the single-particle coordinates 
of the quarks, (r 1; r 2 , . . . , r 5 ), as usual. The angular motion of the system is described 
with a solid spherical harmonics ye m (r) = r e Yi m (r). Here a positive-definite, symmetric 
matrix A and a vector u are parameters which define the shape of the orbital part of the 
basis wave function. Note that u serves to define an appropriate coordinate responsible for 
the rotational motion of the system and that the spherical part of the orbital function can 
equivalently be expressed in terms of the relative distance vectors of the quarks: 



exp 



-—xAx 
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exp 



{ - Y a ki( r k - n) 2 } 



(9) 



l>k=l 



with appropriate a^'s. We use the correlated Gaussian because there are many exam ple s 



22, m 



which demonstrate its power for an accurate description of few-particle systems 
Clearly our orbital part is translation-invariant, so our theory is free from any spurious 
center of mass motion. Note that the parity of the G + is given by P = (-1) L+1 in the 
present formalism. 

One of the most natural coupling schemes for the spin, isospin and color parts is to follow 
a successive coupling. The spin part may be expressed as 



,/,( s P in ) 

rSMs 



sc 

^(512512351234)57^3 



[ki(l)xi(2)]s uX i(3)L v 4 (4) c „ xi(5) 



5l234' 



(10) 



SM S 



where Xim s (0 i s ^ e s Pi n function of the ith quark. Table El lists possible sets of 
(S12, S123, S1234) for S = I and |. For the isospin part we assume T = 0, Mt = 0. Be- 
cause the s does not carry an isospin, the isospin part can be given by coupling four ~ 
angular momenta: 



.(isospin) ^SC 
VTMt ~ S(Ti 2 Ti2 3 0)T=0M 7 



[[ei(i)ei(2)]T 12 ei(3)] Tl23 ^(4)] oo e(5), (n) 



6 



TABLE II: Intermediate spin labels in different models 



9 - I 
° - 2 
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Successive coupling Diquark-diquark type Baryon-Meson type 

gig J 123 <->1234 J 12 ^34 J1234 g!2 ^123 ^45 

i 00 i 
Oil 01 1 Oil 

1 i 11 1 i 

1 \ 1 10 1 1 i 1 
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Successive coupling Diquark-diquark type Baryon-Meson type 

coo o o o o o c 

gig g!23 <->1234 gl2 »->34 J1234 J 12 J 123 '-'45 

1 | 2 11 2 1 | 

1 I 1 1 1 1 I 1 

0^1 01 1 Oil 



2 " "2 

I 1 1 1 1 1 I 



where £i mt (i) is the isospin function of the zth quark and £(5) the "isospin" function of the 

s. Possible values of (Tj^lm) that enable one to make T = 0, Mt = are (0, |) and (1, ~). 
Similarly, the color part can be given, in the successive coupling, as 



./.(color) ^SC 

^(00)000 Ly (ri 2 r 123 r 1 234)(oo)ooo 



[[C (1 o)(l)C (1 o)(2)] ri2 C (1 o)(3)] ri23 C {10) (4)] r C {01) (5) 



, (12) 

(00)000 



where C(i) stands for the color function of the ith quark, and the square bracket [CriCr 2 ]r 
denotes the SU(3) coupling of two functions with SU(3) irreducible representations Ti and 
T2 to that of a definite SU(3) representation T. Here 1^1234 must be (10) to make a color 
singlet + . Possible sets of (Pi2, I123, Ti234 = (10)) are listed in Table ITTT1 As seen from the 
table, three channels make a complete set for the color space of the + . 
To sum up, we may express a trial wave function for the + as 

/C 

*;m = E^. (13) 

-i=i 

with the basis function 

$i = aj)x(s 12 s 12 3Sfi234)s]jAf^(Ti2T m o)ooC(ri 2 ri23ri234)(oo)ooo|, (14) 

where i stands for a set of (A, u, Si 2 , S123, S1234, T 12 , T 12 , r^). Here T123 and Ti234 are 
omitted as they have to be \ and (10), respectively. Note that the 4x4 matrix A is specified 
by ten parameters (An, A 12 , . . . , A 44 ) or equivalently («i2, QJi 3 , . . . , a 45 ) and the u by three 
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TABLE III: Intermediate color SU(3) labels in different models 



Successive 


coupling 


Diquark-diquark type 


Baryon-Meson type 


Ti 2 


Tl23 


Tl234 


Tl2 


Tl234 


Tl2 Ti23 


(01) 


(00) 


(10) 


(01) (01) 


(10) 


(01) (00) (00) 


(01) 


(11) 


(10) 


(01) (20) 


(10) 


(01) (11) (11) 


(20) 


(11) 


(10) 


(20) (01) 


(10) 


(20) (11) (11) 



parameters (as it can be set to a unit vector, uu — 1, without loss of generality). Thus we 
have 13 (10) continuous parameters for the basis function with L^0 (L = 0) and in addition 
the spin-isospin-color channel label which is one of 5 x 2 x 3 = 30 channels for S = \ or 
4x2 x3 = 24 channels for S— |. We stress that all the possibilities satisfying the prescribed 
quantum numbers (S, T = 0, T = (00)) can be taken into account in the present formalism. 
Once the basis functions are set up, the coefficients C{ and the eigenvalue M of the 
Hamiltonian can be determined by solving the following eigenvalue problem: 

1 £[(* t \H\* j )-M{* i \* j ))C j = 0. (15) 

The label SC is omitted from the spin-isospin-color parts in Eq. (|14|) because other 
coupling schemes may be used equally well as explained below. 



Transformation of basis 



The successive coupling scheme introduced in sect. Ill Bl is systematic in its construction. 
The basis functions in that scheme constitute a complete set for specified quantum numbers 
and make it possible to represent any coupling schemes of physical interest. 

The + is considered a system of two diquarks and the s in the diquark model 0, so it 
is useful to define the following coupling scheme, e.g., for the spin part 



MSi2S 3 4S 123 i)SM s 



l\Xh (l)Xi (2)]s 12 \Xi (3)X4 (4)kJs 12 , 4 Xi (5) 



SM S 



(16) 



Here one diquark has spin S 12 and the other diquark S34, and they are coupled to S1234 . 
We list possible sets of (S'12, S34, 61234) in Table EJ As mentioned, this spin function can be 
expressed in terms of combinations of the spin functions in the successive coupling. This 
expansion can be represented as follows (by replacing the spin label S with J): 




^(^12^3^1234^4; ^123^34) 
J123 j x 



JM 




JM 



(17) 



S 



Here a Racah coefficient in unitary form, U, is expressed in terms of the 6j symbol 
U(J 1 J 2 JJ 3 ;J 12 J 23 ) = ((JiJ 2 )Ji2, Js;JM\J u (J 2 J 3 )J 23 ; JM) 

v''-^-'-*\/(2j,, + iu2.7.,, n/ Jl 1/2 Jl2 l. 

■h J J23 J 



(18) 



To be more explicit, the spin function in the diquark model is transformed to the one in the 
successive coupling scheme by the matrix 



/l 00 

1 

1 











J\ 



2 
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1/ 



for S=-, 



(19) 



and 
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1 
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/ 



for S 



3 
2' 



(20) 



Here the spin functions in both models are arranged in order of those defined in Table ITT1 

Similarly, the color function in the diquark model can be expressed in terms of the suc- 
cessive coupling scheme. Possible sets of (r 12 , r 34 , Ti 234 ) in the diquark model are listed in 
Table 11111 The transformation between the diquark coupling and the successive coupling 
can be performed using the formula (1171) (by replacing the Racah coefficient with the corre- 
sponding one in the SU(3) algebra 28, 29]). Here the angular momentum label J should be 
understood to denote the color SU(3) label Y. Note that all the SU(3) couplings appearing 
here are multiplicity-free. The color function in the diquark model is expressed in terms of 
the one in the successive coupling through the matrix 



(- 



\ 



- I- 

3 V 3 

3 V 3 

1/ 







(21) 



A baryon- meson (q 3 -qq) coupling scheme is also important. Here three among the four 
quarks are coupled to a state with a definite spin, isospin and color, and the remaining quark 
and the s are coupled to a state with given quantum numbers, and finally they are coupled 
to a resultant state. Again we can show this coupling scheme in a pictorial way and express 
it in terms of the successive coupling: 




JM 



J5 — X/ ^(^123^4^5) J1234J45 
J1234 




JM 



(22) 
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Possible sets of (J12, J123, ^45) values in the baryon-meson coupling are listed in Table HP for 
the spin part (x BM ) and in Table ITTT1 for the color part (C BM ), respectively. Note that 
and are not necessarily (00), that is, a colored baryon and a colored meson must be 
included to obtain a complete basis for a color singlet pentaquark. The spin function in the 
baryon-meson model is transformed to the one in the successive coupling by 



/ _l yj 

2 

yj 

2 







yj 

2 
1 
2 










_ 1 

2 

yj 

2 








yj 
2 
1 

2 





o\ 







forS = -, 



(23) 



and 



V 













o\ 



1 

1/ 



for 5* 



3 
2' 



(24) 



The color function in the baryon-meson model is transformed to that of the successive 
coupling by a unit matrix. 

The Karliner-Lipkin model proposed in |l0| considers the + as a coupled system of a ud 
diquark and a uds triquark. Here the ud pair in the diquark has Su = 0, Ti 2 = 0, ri 2 = (01), 
and the other pair in the triquark has 5*34 = 1, T34 = 0, 1^34 = (20). The quantum numbers of 
the triquark is S345 = ~, T 345 = 0, r 345 = (10). Recoupling this configuration to the diquark- 
diquark model, we find that the Karliner-Lipkin model is reduced to one particular basis in 
the diquark model with S1234 = 1, T 1234 = 0, r 12 34 = (10), namely 




The calculation of the decay width of the B + will be reduced to an overlap of the + wave 
function with an nK + channel wave function. See sect. UVBl Here we write the spin-isospin- 
color part of the nK + channel wave function in terms of the successive coupling scheme. 
Both n and K + have isospin 1 and only the isospin-singlet component of the coupled state 
contributes to the overlap with the + wave function. The nucleon (N) is symmetric in its 
orbital space as well as in the spin-isospin space. Combining the transformation matrices 
given above, we can express the spin-isospin-color part of n and K + , ip(n)ip(K + ), in terms 
of the successive coupling as follows (by omitting the isospin-triplet component): 

+ ( - 2 X S°)2 m + ^^a C il)|m)^(l C iO)00 W( S (01)(00)(10))(00)000- ( 26 ) 
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Here the first — -^75 factor comes from the Clebsch-Gordan coefficient of the isospin-singlet 
coupling and the second ^ factor from the symmetric spin-isospin function of the nucleon. 

We briefly comment on the transformation of the orbital part. The coordinate appropriate 
for the diquark model is defined as 

y 1 = r 1 -r 2 , y 2 = r3-r 4 , y 3 = -(n + r 2 ) - - (r 3 + r 4 ), 

Va = \{ri + r 2 + r 3 + r 4 ) - r 5 . (27) 

An orbital function in the diquark model may be expressed as c/) LMl (A', u', y). Because y is 
related to a; by a linear transformation, y = Tx, with an appropriate 4x4 matrix T, this 
orbital function can be reduced to the one in the x representation: 

cj>LM L (A', u, y) = (f) LML (A' T , u' T , x), (28) 

with 

A' T = TA'T, u'rp = Tvl. (29) 

Here T is a transposed matrix of T. Therefore, we do not need to introduce a new diquark 
orbital function explicitly but can take into account it by simply adopting A' T and u' T as 
A and u in Eq. (|8*|). It is clear that an orbital function appropriate for the baryon-meson 
model can be expressed in terms of the correlated Gaussian (JBJ) as well. 

It should be stressed that the present formalism can take into account all possible types 
of correlation in a single framework, so to evaluate matrix elements can be reduced to that 
in a particular representation [23]. That is, one does not need to calculate matrix elements 
separately in the different coordinate sets. 

The calculation of matrix elements can be performed separately in the orbital, spin, 
isospin and color parts. The permutation P of the antisymmetrizer A causes a linear trans- 
formation of the coordinate x to x' in the orbital function; x' = Vx. Thus the matrix 
element of an operator O acting in the orbital space can be calculated as 

((j)LM L (A,u, x)\OP\(J) L >m' l (A',u, x)) = {<f> LML {A,u,x)\0\<f> VWL {VA!V,Vu',x)). (30) 

Here the action of the permutation results in just renaming A 1 as VA'V and v! as Vv! , so 
that we do not need to change the functional form of 4>t, m t a t all. The orbital matrix element 
was evaluated using the method explained in detail in [23(. The matrix elements involving 
the spin, isospin and color functions are evaluated using the Wigner-Eckart theorem as well 
as SU(2) and SU(3) recoupling techniques j28j. 



III. RESONANCE IN A SINGLE-CHANNEL CALCULATION 

The mass of the + is in the continuum above the N + K threshold. Here we briefly 
discuss how to predict a resonance mass using the basis functions explained in the previous 
section. There are some powerful methods such as a complex scaling method [30] and an 
analytic continuation in a coupling constant (3]], |32| to use a generalization of bound-state 
problems for locating the resonance. We use the real stabilization method j2jj among others 
for its simplicity. In the stabilization method the Schrodinger equation is solved in a box, 
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i.e., on a square-integrable basis, which makes all solutions look like bound states. From 
among the bound-state-looking discrete states, those corresponding to the resonances are 
singled out by exploiting their stability against changes of the box size (which is practically 
equivalent to the basis dimension in the present case). 

We illustrate this procedure in a single-channel calculation where the basis functions are 
limited to one particular spin-isospin-color configuration. As the single channel we adopt two 
representative ones: One is the Jaffe-Wilczek model in which both diquarks are restricted 
to 5 = 0, T = 0, T = (01) and these identical bosons are coupled to be antisymmetric in the 
color state, and the other is the Karliner-Lipkin model |l0j in which the triquark with 
S — |, T — 0, r = (10) and the diquark with S = 0, T = 0, T = (01) are coupled to a color 
singlet state (see Eq. (J23j0 . In the single-channel calculation, only the matrix A and the 
vector u are variational parameters characterizing the basis function, and the parameter 
space for them are first defined appropriately. Figure plots the mass eigenvalues vs the 
basis dimension for J p = | + (L — 1, S — |) in the Jaffe-Wilczek model, where the basis 
was randomly chosen from the parameter space without any selection procedure and just 
increased one by one. A remarkable feature of this figure is that one has a solution whose 
mass stays nearly constant. The stable solution is not a ground state but has a large mass 
above the N+K threshold. The rms radius calculated from the stable solution is also stable: 
The mass and the rms radius are 2370 MeV and 0.72 fm for the basis dimension /C=100, and 
2340 MeV and 0.73 fm for /C=200, respectively. This stability is an indication required for 
the existence of a resonance. To extend this type of calculations to a full coupled-channels 
problem is hard, however, because the number of channels is fairly large. It is advisable to 
confirm the stability even though the basis dimension is truncated through some selection 
procedure. 

Examples of such a truncated calculation are presented in Fig. El where the mass as well 
as the rms radius are displayed as a function of the basis dimension. Here the basis selection 
was performed using the SVM [i^, as it leads to virtually exact solutions in rather 
small dimension and its power is proved in diverse few-body problems. In the SVM the 
basis dimension is increased one by one by testing a number of candidates which are chosen 
randomly. The importance of each candidate is evaluated by calculating the energy gain 




12 




FIG. 2: The mass and rms radius of the + in a single-channel calculation. The figures on the left 
side are the results in the Jaffe-Wilczek model, while those on the right side are obtained in the 
Karliner-Lipkin model. The AL1 potential is used. 

produced when it is included in the basis set. The candidate which produces the lowest of 
the energies is chosen to be a member of the basis set. Furthermore, it is shown in |33| that 
the basis set for describing a resonance can be selected very efficiently with the SVM. The 
basis selection in the present calculation was actually done as follows: each element of A and 
u was randomly generated twenty times, and the best one among these trials was selected as 
a tentative candidate for a new basis. This process was repeated ten times to select a most 
successful one from among these tentative candidates. A typical computer time needed to 
determine twenty basis states was about half an hour on a computer loading the Intel (R) 
Pentium (R) 4 CPU 3.20 GHz. 

We see two characteristics in either model of the figure. One is that the lowest mass 
decreases as the basis dimension increases, reaching even a negative value. The other is the 
appearance of such a unique mass eigenvalue that is stable against the increase of the basis 
dimension. By looking at the figure of the rms radius, we note that the radius corresponding 
to the former solution is rather large and continues to increase with increasing dimension. 
In contrast to this case, the radius for the latter solution stays approximately constant, and 
it decreases a little from the one of Fig. ^ According to the criterion of the real stabilization 
method, we may identify the latter mass as the resonance mass and the corresponding wave 
function as the approximate resonance wave function which is valid except at large distances. 
Of course the mass obtained here is just a result of the single-channel calculation, so before 
declaring it to be the + mass we have to consider the effect of other channels which may 
contribute to the resonance. In particular, the effect of the NK channel must be examined 
carefully. This will be performed in the next section. 
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We have performed a single-channel calculation for the other configurations as well in 
each case of the successive coupling, the diquark model or the baryon-meson model. The 
falloff of the mass and the increase of the rms radius are observed in all the cases. The reason 
for this seems due to the color dependence of the quark-quark potential as was pointed out 
in sect. Ill Al First we notice that the mass of the colored subsystem happens to be negative 
and its radius gets very large obviously because of the (\f ■ A?)ry operator which works to 
deconfine the system in special color channels. A five-quark system usually contains such 
colored components, so the variational ground state will make use of it to gain the energy. 

It is nevertheless fortunate that we have the sign of a resonance, that is, the stability in 
the mass diagram, which is indispensable for getting a resonance in the bound-state-looking 
calculations. Those cases which show the stability are characterized by that at least one 
of the diquarks is antisymmetric with respect to the simultaneous interchange of the spin, 
isospin and color degrees of freedom. 

A single-channel calculation has also been performed for other J p states. We have ob- 
tained some stable masses around 2500, 2700 and 2900 MeV for (L = 1, S = §), (L = 2, S = |) 
and (L = 2, S= |) cases, respectively. It seems that the larger L and the larger S the system 
has the larger its mass is. This can be expected from the role played by the kinetic energy 
and the color magnetic piece of the quark-quark potential. We will focus our attention on 
the J P = \ and | states in what follows. 



IV. RESULTS 

A. Mass spectrum 

We have made a coupled-channels calculation to predict the mass of the + for J p = 
f(L = 0,S = §,T = 0), I + (| + )(L = 1,S = i,T = 0) and §~ (L = 0, S = §, T = 0). As 
discussed in sect. Ill Bl there are in principle 30 channels for S = |, T = and 24 channels 
for S = |, T = 0. Among these some channels play an important role to produce a resonance 
but some others play much less significant roles. In the actual calculation we first single out 
basis functions from those channels which give a stable mass in a single-channel calculation, 
and include them in the coupled-channels calculation. Included are 7, 5 and 12 channels for 
J P = \ , | and | + (| + ), respectively. These channels are expressed using the channel labels 
of S12, Tu, r 12 , 5*34, T 34 , r 34 , S1234 in the diquark model as follows: 

for J p = i~ 00(01)00(01)0, 00(01)10(20)1, 01(20)11(01)1, 11(01)11(01)0, 
10(20)00(01)1, 11(01)01(20)1, 11(01)11(01)1, 

for J p = |" 10(20)00(01)1, 11(01)01(20)1, 00(01)10(20)1, 01(20)11(01)1, 
11(01)11(01)1, 

for J p = \ + 00(01)00(01)0, 00(01)10(20)1, 01(20)11(01)1, 11(01)11(01)0, 
10(20)00(01)1, 11(01)01(20)1, 10(01)10(20)1, 10(20)10(01)1, 
11(01)11(01)1, 11(20)11(01)1, 10(01)00(01)1, 00(01)10(01)1. (31) 

To obtain a final result for a resonance, we further consider the effect of those basis 
functions which are generated from the NK channel for J p = | or from the NK* channel 

for I because the coupling of the resonance with that channel appears to be important for 
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calculating the width. Note that the NK channel does not couple to | states as they differ 
in the total spin S, but the NK* channel can have the same spin S = |. For example, the 
NK channel basis functions included are expressed as 

MNK)=a{[[* h (N)%.(K)^^ (32) 

where ^> I= i T= i(N) and ^ I=QT= i(K) are respectively the wave functions of N and K which 
are obtained by solving the three-quark and quark-antiquark Hamiltonians with the potential 
of Eq. (J2J) or Eq. (jSJ). These are coupled to the relative motion function of Gaussian form with 
the orbital angular momentum I to the total angular momentum JM. Here the coordinate 
z 4 = Rn — Rk is the relative distance vector between the center of mass of N, R^, and 
the center of mass of K, Rk- In the calculation the length parameter is set to be \j J~Oi — 
0.6xl.4 i - 1 fm(i = l,2,...,5). 

The Hamiltonian is diagonalized in the set of the chosen basis functions. Figures El 
and |U display the calculated mass as a function of the basis dimension. We clearly see a 
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FIG. 3: The mass of the + in a coupled-channels calculation. The NK channel is explicitly 
included for J p = \ , while the NK* channel is for | . The AL1 potential is used. 
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FIG. 4: The same as Fig. El but for the TS potential. 



stabilized mass in the diagram and have confirmed that the rms radius calculated from the 
corresponding wave function is stable. This stabilized solution is accepted as a resonance we 
seek. We have found that the mass and the radius of the stabilized solution remains basically 
unchanged by including the NK (and NK*) channel of Eq. (j32j) . Evidently, our resonance 
does not correspond to a variational energy minimum. We have used a variational method 
in order to set up compact, appropriate basis functions. This is in contrast to the approach 



m 



111 . 12l | . where a variational energy minimum is looked for and moreover a coupling to 



the baryon-meson channel like NK is excluded from the model space. 

The masses and rms radii obtained in the present calculation are summarized in Fig.|S]and 
Table HVl It is concluded that the mass of the G + increases in order of J p 



1- 3- 1+/3- 

2 '2 '2 ^2 h 

state than for 



which agrees with the result of Ref . |ljj . A lower mass is obtained for the | 

the | + state, which is consistent with the recent variational Monte carlo calculation [13] but 
is opposite to the level order obtained in 0|. The color magnetic and kinetic energy terms 
in the Hamiltonian play the most important contribution to generate the mass differences. 



The mass splitting between the § and 



states is smaller in the AL1 potential than 



that in the TS potential, but both potentials otherwise give qualitatively similar results. Our 
result for the TS potential is compared to that of 11], and we see that the mass splitting is 
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FIG. 5: The mass spectrum of the + from the N+K threshold. The result of Ref. ll| is due to 
the TS potential. 



not necessarily the same between the two. This is probably due to the difference in obtaining 
the mass as well as in specifying the correlation of the quark dynamics. 

The calculated mass of the B + is higher than the observed mass (~1540 MeV) by about 
500-600 MeV. This appears to be common in all the calculations using the potential of 
OGE type 11, 0, 3 o r the OGE plus one-pion exchange potential ^3(. The zero-point 
energy term is uncertain, so it is hard to get an absolute mass unambiguously. There are 
some discussions to reduce the mass. One is to consider the flavor-spin dependence of the 
quark-quark potential or the one-boson exchange potentia l |l5l Il7l| . The second is to take 
into account the role of the instanton induced interaction |3J,|35J]. Though this interaction 
is expected to bring about an extra attraction, its effect has not yet been quantified. The 
third is to consider the relativistic effect of the quark motion. An extensive calculation has 
been made in 
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which used the semirelativistic expression for the kinetic energy and the 
one-boson exchange potential as well as the OGE potential. It seems that the | + state 

comes down close to or even lower than the § state. The calculated mass is, however, still 
larger than the observed mass, and even goes around 2 GeV if the zero-point energy is taken 



to be proportional to the number of quarks. 



TABLE IV: The mass and rms radius of the + . The calculated N+K threshold energy is 1486 
MeV for the ALl potential and 1451 MeV for the TS potential, respectively. 
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B. Decay width and spin-parity 



The width of a resonance is in principle calculated through a coupling to the continuum 
states of decay channels. As we have localized the 6 + , it may be appealing to use, for 
example, the complex scaling method to calculate the width as was done in 33[. This is, 
however, too sophisticated in the present case since the calculated 6 + mass is subject to 
large uncertainty of the quark-quark potential. On the other hand, the wave function of 
the resonance remains the same for an arbitrary adjustment of the resonance energy due to 
the change of the zero-point energy term, so we instead use the .R-matrix theory [3a . |3?| i n 
which the width T of an isolated resonance is calculated through 



T = 2P,(a) 7 2 (a), 

where a is a channel radius and Pi{a) is the penetrability given by 

PA a) 



ka 



ilika) + nf{ka) 



(33) 



(34) 



Here i and k are the orbital angular momentum and the wave number of the relative motion 
between the decaying particles, and ji and ri£ are the spherical Bessel functions. When both 
of the decaying particles are charged, the spherical Bessel functions should be replaced by 
the Coulomb functions. 

Most crucial in Eq. (}33|) is the reduced width 7 2 (a), which is related to the reduced width 
amplitude y(a) for the decay to a baryon B and a meson M: 



7 W = —y [a) 



y[r) 




IbT b (B) ^i mTm {M)]i tm t Y e (z 4 



5(z4 — r) 



JM 



z 4 r 



JM 



(35) 



Here /i is the reduced mass of the baryon and the meson, ^>i b t b {B) is the baryon wave 
function which is properly antisymmetrized and normalized, and ^i m t m (M) is the normal- 
ized meson wave function. The coordinate z 4 is the relative distance vector between the 
baryon's center of mass and the meson's center of mass. Introducing a complete basis set 
{/ n ^(r)Y£ m (f)}, the calculation of y{r) is reduced to that of the overlap: 



4' /r 

Try E fne(r){[[^i B T B (B)^ lMT AM)}iTM T fne(z4)Y £ (z^} JM 

n 



JM, 



(36) 



Figure El displays the reduced width amplitude of the + resonance obtained with the 
AL1 potential. The spins of the decay channel are taken as = (§,0) for J p = | and 



-2 - 



1) for J 1 



2 ' 



respectively. Note that the + with J p = | cannot decay 



to the nK + channel because it has S = | and its reduced width amplitude to that channel 
vanishes. For it to decay to the nK + channel, one has to take into account tensor components 
of the interaction between the quarks. However, the | state can decay to the nK* + channel, 
so the reduced width amplitude for this decay with (J, £) = (|, 0) is displayed in the figure. 
It is found that the reduced width amplitude hardly changes with the inclusion of the NK 
(or NK*) channel basis functions defined in Eq. (J32j) . The amplitude of the Wigner limit 
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(WL), y(r) = 




is also plotted in the figure. The Wigner limit is the partial width of a 



particular state whose wave function is uniform up to r and zero beyond, and it is a measure 
to judge whether the resonance has a large component in the decay channel or not. Since 
the calculated reduced width amplitudes are considerably small compared to the Wigner 
limit, we can conclude that the + does not have large nK + component. This point will be 
confirmed in sect. IIVCI It is noted that the amplitude of the | state is especially small at 
r > 1 fm. 

We estimate the following decay width according to Eq. (}33|) : the nK + decay of the G + 
with ~ and the nK* + decay of the 9 + with | . Because the calculated mass is subject to 
change due to the uncertainty of precise knowledge on the interaction between the quarks, 
we display in Fig. [7| the width as a function of the decay energy E. The channel radius a 
relevant to the decay is estimated as a sum of their radii, and it is about 1.1 fm using the 
values of Table HI The channel radius dependence is also presented around a = 1 fm. 

First let us accept the calculated mass as it is. Then we conclude that no pentaquark 
states G + appear around 1540 MeV, but their masses are expected to be larger than 2 GeV. 
The decay width is then very large, which is different from the conclusion of [l4j. Even 
for the | + state, the width would be in the order of 100 MeV. In the case of J p = | , the 
mass from the n+K* + threshold is 270 MeV for the AL1 potential and 320 MeV for the TS 
potential, respectively, so that the decay width to this channel is a few MeV. The magnitude 
of this width is probably small enough to be observed as a sharp resonance. 

Next let us assume that, by shifting the calculated mass spectrum, one of the calculated 
@ + resonances corresponds to the observed mass, 1540 MeV, namely 100 MeV above the 
n + K + threshold. There are following three possibilities, and we estimate the decay width 
from Fig. [7| for each case: 




0.5 1 1.5 2 2.5 

r [fm] 

FIG. 6: The reduced width amplitude of the + to the nK + channel for J p = ^ and to the 
nK* + channel for J p = | . The reduced width amplitude for J p = | + is the same as that for 
J p = 1 + . The AL1 potential is used. 
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(1) If the spin-parity of the G + is ~ , its width is several tens to 100 MeV, which is too 

large to be compared to the observation. 

(2) // the spin-parity of the G + is | , its nK + decay width is practically zero. Because of 

the mass spectrum predicted by the present model (see Fig El), another 6 + with \ 
is expected to appear below the + (§ ) as well and its width is in the order of a few 
MeV. 

(3) // the spin-parity of the + is | + (§ + ), its nK + decay width is about 10 MeV. Below 

this resonance other 6 + states with | and | appear as quasi-bound states. 

The case (1) is apparently in contradiction to experiment. The case (2) may remain as an 
actual possibility because the decay width is practically zero. However, in this case another 
pentaquark state with J p = | with a small width should also be observed near or slightly 
above the n+K + threshold. No such experimental information is available at present. The 
case (3) can also be an actual possibility and the decay width is in the order of 10 MeV. 
Though the existence of the 9 + states with \ and | are also predicted in this case, their 




100 200 300 400 500 600 700 100 200 300 400 500 600 700 

E [MeV] E [MeV] 




FIG. 7: The decay width V of the + to the nK + channel for J p = | and to the nK* + channel 
for J p = | . E is the decay energy from the respective threshold and a is the channel radius. The 
AL1 potential is used. 
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energies are below the n + K + threshold. Thus we expect no signal for the existence of 
such states from the type of experiment which uses the K + n invariant mass spectrum for 
identifying the pentaquark. Therefore, the case (3) that the spin-parity of the G + (1540) is 
| (| ) seems not to be in contradiction with the available experimental information. 

C. Structure 

To discuss the structure of the 6 + , it is useful to calculate a single-particle density and 
a correlation function between the particles. These functions F are all defined as 

F{r) = ^ p JM \8{ux-r) \^ P JM ), (37) 

where u = (ui, uj 2 , w 3 , U4) is a vector which is chosen appropriately depending on the quan- 
tity of interest. For example, the density distribution of the ud quarks from the center of 
mass -R c .m. of the system, the quark-quark correlation function and the quark- ant iquark (s) 
correlation function can be calculated by choosing u— (0, 0, — |, im m + m ) (&x = r i —R c . m ) 1 
uj — (1, 0, 0, 0) (u}x = ri — r 2 ), andu; = (0, 0, — |, 1) (ux = r 4: — r 5 ), respectively, because the 
wave function is properly antisymmetrized. The choice of uo = (0, 0, 0, 1) gives the distribu- 
tion (q 4 -s) of the s from the center of mass of the four quarks. When the delta function, 
5{ojx — r), is expanded in multipoles, only the monopole term contributes to the function 
F(r) for both cases of J p = i , so F(r) becomes spherically symmetric. 

Figure |H1 displays the density distributions as well as the correlation functions for the G + 
with J p = I . It is seen that the distribution of the s is confined in the smaller region around 
the center of mass of the system than the uci-quark distribution. The probability density 
reaches a maximum at about 0.30 fm for s and about 0.44 fm for ud, respectively. This 
feature seems to arise from the fact that the quark-antiquark interaction is stronger than the 
quark-quark interaction. Corresponding to this feature, we see from the correlation function 
that the average distance between the ud and s quarks is shorter than that between the ud 
quarks. Figure H3 displays the similar distributions for the + with | + . The distribution 
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FIG. 8: The density distributions of the quarks (left) and the correlation functions between the 
quarks (right) for the + with J p = ^ . The AL1 potential is used. 
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FIG. 9: The same as Fig.EJbut for J p = \ + . 

extends to larger distances than that of J p = ~ . Again, the s distribution has a peak at 
smaller distances from the center of mass than the ud distribution. 

Another structure information is obtained by decomposing the wave function into various 
channels of the spin, isospin and color spaces. To this end, let us define the channel wave 
function $ c . See sects. HTTJ1 and ITTO In the diquark-diquark decomposition $ c is defined as 

<f> _ 1 J .J)D <rDD ^yDD 

- f— : ~S A(Sl2S , 34S , l234)S , M S ^(Ti2T340)00 O (ri2r 3 4(10))(00)000 

^2(1 + d 12 , 3 4> I 

A(5 34 s'i25i234)5M s ?(T34ri20)oo L '(r34ri2(io))(oo)ooo h 
where #12,34 = ($Si 2 S3 4 ( $Ti2T34( s >ri2r 3 4, and cr is the phase defined by 

q _ I ^•Si2+5'34 — 5l234+7 1 i2+T34+Ai2+A34+/^i2+^ 3 4 — 1 (39) 

The channel index c is characterized by a set of values (Si 2 , T 12 , T 12 , S345 T 34 , r 34 , S1234, n), 
where 7r = ±l determines the parity of the $ c . The $ c with 7r = +l is symmetric with respect 
to the simultaneous interchange of the spin, isospin and color parts of the two diquarks, and 
the $ c with 7r = — 1 is antisymmetric. In the case of ($12,34 = 1, a combination of na = 1 only 
is possible, that is, either 7r = +l or it = — 1 is allowed depending on whether Si 234 = 1 or 
5*1234 = 0. In the baryon-meson decomposition $ c is defined as 

WcMs - A(Si25i2 3 545)S , M s ^(Ti2Ti23|)00 ( - / (ri2ri23r 4 5)(00)000- l 4U J 

Here the channel index c stands for (Si 2 , T 12 , T u , S123, T i23 , r i23 , 6*45, r 45 ). In particular, the 
NK channel wave function is defined as 

v N km s - -^|X( lo)iM s ^(Oil)00 + X(ii )iM s ?(iii)oo|°((oi)(oo)(oo))(oo)ooo- l 4i J 

The probability P C (J P ) of finding the channel c in the + is calculated as the expectation 
value of the projector J2m s I ®cM s )(®cM s l : 

P C (J P ) = J2(*JM I ®cM S ){®cM s I *^r>. (42) 
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Values of P c {\ ) are listed in Table IVI for the diquark model and in Table IVT1 for the 
baryon-meson model. The channel of the first row in Table IVI represents the diquark model 
of Jaffe-Wilczek [9f. It is seen that this channel occupies a relatively large component in the 
| state as expected but its magnitude is not overwhelmingly large. Following this channel, 
the diquark pair with SnT^ 12 = 11(01) and 534X341^34 = 11(20) occupies a significant weight 
for J F = i + . The channel of the third row corresponds to that of the diquark-triquark 
model as shown in Eq. ([25|) . We see that the configuration of the Karliner-Lipkin model 
is dominating in the | state but it is not a main configuration in the ~ state. With 

the probability of more than 80 %, the two diquarks in the \ state are found in either the 
Karliner-Lipkin channel or the channel with Si 2 T 12 T 12 = 11(01), S^T^T^ = 11(01), S1234 = 1. 
It is seen from Table IVD that the components in the baryon-meson decomposition are spread 
over many channels. Relatively large components do not necessarily appear in colorless 
baryons. The probability summed over the colored baryon-meson channels is 67 %, which 
signals the importance of the hidden color components. The NK component, Pnk(J P ), is 
calculated to be 0.14 for the | state and 0.073 for the \ + state, respectively. The magnitude 
of this channel is related to the reduced width to the nK + channel discussed in sect. IIVB1 
and it is certainly small, especially for the | + state. 

V. SUMMARY 

We have studied the mass and the decay width of the + consisting of uudds quarks in 
the constituent quark model. The quarks are assumed to interact via a variant of the one- 
gluon exchange potential plus a phenomenological confinement potential which reproduce 
the masses of e.g., N, K and K* reasonably well. The five-quark states are described using 
the correlated Gaussian basis together with the general spin-isospin-color wave functions. 
The basis functions are set up with the stochastic variational method. One of the advantages 
of the present approach is that the importance of all the possible color configurations can 
be tested, so that the hidden color components are naturally taken into account in the 
calculation. We have performed coupled-channels calculations to predict the + masses for 
J p = f(L = 0,S= |,T = 0), I + (f + )(L = l,»S=±,T = 0) and §" (L = 0, S = §, T = 0). 
The NK and NK* channels are explicitly included as they are expected to be important 
to obtain the decay width. The real stabilization method is successfully used to locate the 
resonance in the continuum states above the N+K threshold. The ground state is found to 
be \ and its mass is about 2 GeV, which is 500-600 MeV higher than the observed value. 

The first and second excited states are § and | + (| + )- The color magnetic and kinetic 
energy terms of the Hamiltonian are the most important pieces which contribute to this 
level ordering. 

To estimate the decay width of the G + , we have used the i?-matrix theory in which 
most crucial is the reduced width amplitude for the relevant decay channel. The decay 
channel is specified by the angular momentum /, the addition of the spins of the baryon 
and the meson, as well as the relative orbital angular momentum i between them. The 
calculated decay width is the nK + decay with {1,1) = (|,0) for ~ and (I,£) = (|, 1) for 

| + (| + ), respectively. For the Q + with J F ' = \ , the decay to the nK + channel is forbidden 
because the corresponding reduced width amplitude vanishes, so the decay width to the 
nK* + channel is calculated. All the reduced width amplitudes are small compared to the 
Wigner limit value. As the mass of the + is very large compared to the threshold, the decay 
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width becomes very broad except for the | state whose decay width to the nK* + channel 
is only a few MeV. We have also estimated the decay width by shifting the calculated mass 
to the observed mass, i.e., 100 MeV above the n+K threshold. This shift causes no change 
in the resonance wave functions. The possibility that the observed state is | is ruled out 

because its width is still too large. The case that it is either | or ^ + (| + ) appears not 
to be inconsistent with the observation of the small width. However, in both cases other 
pentaquark state with different J p is expected to exist below the observed one at 1540 MeV. 
Its experimental confirmation will give us a support for the + . 

The structure of the | states is discussed through the density distributions of ud and s 
quarks as well as the two-particle correlation functions of the quarks. It is found that the 
s quark has a narrower distribution near the center of mass of the uudd quarks than the 
ud quark. The average q-q distance is longer compared to that of q-s. These characteristics 
can be understood from the fact that the q-q interaction is stronger than the q-q interac- 
tion. We have decomposed the resonance wave function into various components using the 
diquark-diquark model and the baryon-meson model. The baryon-meson configuration of 
the Karliner-Lipkin model is in fact a main component in the | state but not an important 
component in the | + state, differing from its expectation. In contrast to this, the diquark- 
diquark configuration of the Jaffe-Wilczek model occupies the largest component in the | + 
state as expected, but its magnitude is not very large. The components are actually spread 
over many diquark-diquark channels. 
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TABLE V: Decompositions of the + resonances into the spin, isospin and color channels of the 
diquark model. See Tables [H] and IIIII for the channel labels. T1234 = and 1^1234 = (10) are 
abbreviated, tt is the parity of the channel wave function with respect to the diquark-diquark 
exchange. The AL1 potential is used. 



S12T12 


Tl2 




51234 


TT 


Pc{\ ) 










(01) 








(01) 





-1 


0.0005 


0.2067 








(01) 


1 





(20) 


1 


-1 


0.0017 


0.0457 








(01) 


1 





(20) 


1 


+1 


0.4933 


0.0272 





1 


(01) 


1 


1 


(01) 


1 


-1 


0.0008 


0.0316 





1 


(01) 


1 


1 


(01) 


1 


+1 


0.0008 


0.0316 





1 


(20) 


1 


1 


(01) 


1 


-1 


0.0012 


0.0310 





1 


(20) 


1 


1 


(01) 


1 


+1 


0.1649 


0.0126 


1 


1 


(01) 


1 


1 


(01) 





-1 


0.0003 


0.0062 


1 


1 


(01) 


1 


1 


(01) 


1 


+1 


0.3290 


0.0215 








(01) 








(20) 





-1 


0.0002 


0.0149 








(01) 








(20) 





+1 


0.0002 


0.0149 








(01) 


1 





(01) 


1 


-1 


0.0007 


0.0150 








(01) 


1 





(01) 


1 


+1 


0.0007 


0.0150 





1 


(01) 


1 


1 


(20) 


1 


-1 


0.0000 


0.0123 





1 


(01) 


1 


1 


(20) 


1 


+1 


0.0006 


0.0068 





1 


(01) 





1 


(20) 





-1 


0.0001 


0.0415 





1 


(01) 





1 


(20) 





+1 


0.0001 


0.0415 


1 


1 


(01) 


1 


1 


(20) 





-1 


0.0003 


0.1151 


1 


1 


(01) 


1 


1 


(20) 





+1 


0.0003 


0.1151 


1 





(01) 


1 





(20) 





-1 


0.0001 


0.0415 


1 





(01) 


1 





(20) 





+1 


0.0001 


0.0415 


1 


1 


(01) 


1 


1 


(20) 


1 


-1 


0.0008 


0.0166 


1 


1 


(01) 


1 


1 


(20) 


1 


+1 


0.0008 


0.0166 


1 





(01) 


1 





(20) 


1 


-1 


0.0006 


0.0299 


1 





(01) 


1 





(20) 


1 


+1 


0.0006 


0.0299 





1 


(01) 





1 


(01) 





-1 


0.0000 


0.0000 








(20) 


1 





(01) 


1 


-1 


0.0000 


0.0041 








(20) 


1 





(01) 


1 


+1 


0.0005 


0.0080 


1 





(01) 


1 





(01) 





-1 


0.0000 


0.0000 


1 





(01) 


1 





(01) 


1 


+1 


0.0006 


0.0057 



1.0000 1.0000 



26 



TABLE VI: Decompositions of the + resonances into the spin, isospin and color channels of the 
baryon-meson model. See Tables ITT1 and ITTT1 for the channel labels. The AL1 potential is used. 
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